%================================================
% INITILIALIZE STRUCTURED MESH
%================================================


function [g,p,b,t]=C_unstructured_mesh(domain);



%============================================
% Domain
%============================================
x0=domain(1,1);
x1=domain(1,2);
y0=domain(2,1);
y1=domain(2,2);

%============================================
% geometry
%============================================
g=[
     2     2     2     2
    x0    x1    x1    x0
    x1    x1    x0    x0
    y1    y1    y0    y0
    y1    y0    y0    y1
     0     0     0     0
     1     1     1     1
 ];
 
%============================================
% points
%============================================
p=[
     x0     x1     x1    x0    (x0+x1)./0.5
     y1     y1     y0    y0    (y0+y1)./0.5
 ];


%============================================
% boundary
%============================================
 b=[
     1     2     3     4
     2     3     4     1
     0     0     0     0
     1     1     1     1
     1     2     3     4
     0     0     0     0
     1     1     1     1
 ];
 
 
%============================================
% TRIANGLES
%============================================

 t =[

     2     3     4     1
     1     2     3     4
     5     5     5     5
     1     1     1     1
 ];
